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Abstract. The properties of LRS class II perfect fluid space-times are analyzed using 
the description of geometries in terms of the Riemann tensor and a finite number of 
its covariant derivatives. In this manner it is straightforward to obtain the plane and 
hyperbolic analogues to the spherical symmetric case. For spherically symmetric static 
models the set of equations is reduced to the Tolman-Oppenheimer-Volkoff equation 
only. Some new non-stationary and inhomogeneous solutions with shear, expansion, 
and acceleration of the fluid are presented. Among these are a class of temporally 
self-similar solutions with equation of state given by p = (7 — l)p, 1 < 7 < 2 , and 
a class of solutions characterized by a = — We give an example of a geometry 
where the Riemann tensor and the Ricci rotation coefficients are not sufficient to give 
a complete description of the geometry. Using an extension of the method, we And the 
full metric in terms of curvature quantities. 
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1. Introduction 


It is well known that manifolds with a metric strnctnre, i.e., geometries, are locally 
completely determined by a set consisting of the Riemann tensor and its covariant 
derivatives with respect to a frame with constant metric components H- This set 
will henceforth be called 


■ {R-ijkh Rijkl\mn -} 


where p is such that the components in are functionally dependent on those in 

(as functions of both the coordinates of the manifold and the parameters of the 
generalized rotation group). This description can be used to classify geometries |||,^,|| 
and also to construct new solutions to Einstein’s equations in terms of the set R^^^ |^|-|^. 
A set RP~^^ describes a geometry iff certain integrability conditions are satisfied ||^,|^. 
These are equivalent to the Ricci equations up to order p + 2 and also part of the 
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Bianchi equations if the geometry has isometries (note that the Ricci equations in our 
terminology also contain part of the commntators). If the geometry has isometries one 
does not get the full line-element with this method. In ii it was shown how the 
method can be extended to give the basis 1 -forms uj\ and hence the fnll line-element in 
terms of canonical coordinates, by using the isometry algebra, which is easily obtained 
from |[T0|| . 

In an earlier paper § one of the anthors nsed this description to constrnct locally 
rotationally symmetric (LRS) perfect fluid solutions. The LRS classes I (stationary 
solntions with vorticity) and III (space-homogeneons with spatial twist) where stndied. 
The full line-element could be constructed from in all cases. In this paper we study 
the remaining class II containing solntions that in general depend on both a time-like 
and a space-like coordinate, but have no vorticity or spatial twist. There are subclasses 
of static solutions and spatially homogeneous solutions in this class. The class have of 
course been extensively stndied and is perhaps the most physically interesting 

among the LRS models since it contains, e.g., solutions describing the gravitational 
collapse of protostars (as far as energy dissipation may be neglected) and relativistic 
spherical stars in equilibrium (if once again energy transfer is neglected). In earlier 
work one has often assnmed spherical symmetry. It is worth pointing ont that to many 
snch solntions there are corresponding solntions with hyperbolic symmetry (on the two- 
dimensional surfaces of symmetry). 

The set can be obtained from the smaller set S = {R’‘jki,^%i,x°‘\i}, where 

T'fej {% = l,..,n{n — l)/2— dim(isotropy group)) are the part of the Ricci rotation 
coefficients that generate the rotations that do not leave invariant, (a = 1 ,..., n— 
dim(orbits)) are essential coordinates on the manifold (i.e., the elements in R^ are 
fnnctions only of x°‘ when the frame is fixed) and the coordinate gradients are given 
by x°‘\i := ei{x°‘) where is the dnal basis to uj\ Onr method have similarities 
with the covariant approach by van Elst and Ellis [^], but they use the smaller set 
S' = {R'‘jki,T%i}. In section 3 we give an example of two geometries with the same 
set S', but with different x°]i, resulting in different second covariant derivatives of the 
Riemann tensor. Hence the geometries are ineqnivalent and the set S' is not snfflcient 
in general to give a complete local description of the geometry. In this section we also 
discuss the general structure of the equations and the possible choices of coordinates. In 
section 4 the static case is treated and the set of eqnations is (given an eqnation of state 
R = Rip)) reduced to one integro-differential equation for the pressure gradient (the 
Tolman-Oppenheimer-Volkoff eqnation). Spatially homogeneons models are stndied 
in section 5. In general the system is reduced to two differential equations. Tilted 
models are obtained when the Gaussian cnrvatnre K of the 2-snrfaces of symmetry 
{ dt = dr = 0} eqnals zero. In section 6 , models depending on both a time-like and a 
space-like coordinate are treated. Some probably new solutions with shear, expansion 
and acceleration are given. One class is characterized by a = —(where a is shear 
and 0 expansion). A class of temporally self-similar solutions with iL = 0 and equation 
of state given by p = (7 — l)/i, where 1 < 7 < 2, is also presented. Einally, in section 7 
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2. Preliminaries 

Perfect flnids with locally rotationally symmetry (LRS) may be divided into three 
separate classes |||,|^,|^. In [|] one of us studied the classes I and III. Here we treat the 
remaining class II. The vorticity of the flnid, the spatial twist of the preferred spatial 
direction and the magnetic part of the Weyl tensor all eqnal zero in this class. 

The energy-momentum tensor for a perfect fluid is given by T = {fi+p)uC)u — pg , 
where u is the 4-velocity of the flnid, p the energy density, p the pressnre and g the 
metric. We choose a comoving Lorentz-tetrad, i.e. we select the tangent basis {cj} and 
its dnal {cj*}, i = 0, ...,3, snch that g = rjijOj'^ 0 

where = u. The axes are rotated so that the LRS- symmetry lies in the 23-plane, 
which makes our tetrad fixed up to rotations in that plane. 

As stated above, can be obtained from the set S = Since 

the geometry is assumed to be invariant under rotations in the 23-plane, the T\i are 
exclnded from S. For the same reason there are at most two fnnctionally independent 
elements in RP^^ and all elements can hence be seen as functions of two coordinates, 
say t and r where t is time-like and r is space-like. Their exterior derivatives can be 
written 

Rt = Xuj^ + Yu)^, 
dr = xuj^ + ^ 

so that the tangent basis vectors Cq and ei become 
Co = Xdt + xdr , 

Cl = Ydt + ydr , 

from which the coordinate gradients x°'\i easily can be read off. The reqnirement for the 
coordinates to be independent is Xy — Yxfl^Q. t is time-like when > T^, and r is 
space-like when y'^ < x^. 

The Riemann tensor for a LRS perfect fluid with vanishing magnetic part of the 
Weyl tensor is given by 

-Roioi = E — |(/r -|- 3p) , Ro 202 = Ro303 = ~\E — |(/i -|- 3p) , 

Ri212 = Ri313 = \E — i/r , R2323 = ~E — |/r , 

where E = En and Eif=Cikjiu^u^ = Coioj denotes the electric part of the Weyl tensor. 

Some of the rotation-coefficients are expressible in the kinematic qnantities 
acceleration a*, expansion 0, vorticity uJij (which is zero in our case) and shear (7^ 
through 
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UJij. hi hj ^ij' {l^k) “t“ 7 

Q:='ViU\ ai:=u^'V jUi, 

and hij\=UiUj — rjij is the projection tensor. The choice of a comoving frame results in 
flo = o'io = 0, because of their orthogonality to the fluid 4-velocity. From the requirement 
of LRS there is only one component of the acceleration, ai\=a. The shear will, because 
it is traceless, look like {aifl = diag(0, —2 ct, a, a). The necessary non-zero rotation 
coefficients for LRS class II are then 


Toio — —O' 1 Ton — 2cr -|- |0 a, 

ro22 = Loss = — cr -|- |0 := /3, ri22 = ri33:= — K , 

(see 1^ for the general LRS case). The rotation coefficient k corresponds to the spatial 
divergence of the vector-field Ci. Thus all rotation coefficients are covariantly defined 
because of the LRS, and we have the following set of functions determining the geometry 
in the generic case: 


S = {E, fi, p, a, a, fl, K, X, F, x, y, rjifl . 

We would like to stress the fact that the variables we have chosen to describe 
the geometry are only singled out because of their direct correspondence with the 
components of the Riemann tensor and the Ricci rotation coefficients. Any other set 
that can reproduce the set could be used as well. 

The functions in the set S are subject to the following integrability conditions: 


„/3 _ T^m 

^ |[/c,|/3|^ |/] ^ |m-L [kl] •) 

( 2 , 1 ) 

hij b[j,\a.\^ |i] “T yL “r bk^ [ij] 

. ( 2 . 2 ) 

R [ijk] 0 , 

(2,3) 

R^qlij;k] = 0 . 

(2.4) 

1 , 2 ,...,/, {“b} = 1 , 2 , ...,m, t = / -k 1 ,/ -k 2 , ...,n. 

and {^q} = m -k 1 , m -k 


2,..., n(n — l)/2, where I := n — dim(orbits) and m := n{n — l)/2 — dim(isotropy group) 
(see Refs. ||^,Q for a more detailed exposition of the method). For the case at hand they 
are: 

Commutator equations 

[ei, eo] t = eflX) - eo(F) = XF + X'y -YX - Y'x = Xa + Ya , (2.5) 
[ei, eo] r = eflx) — eo(y) = xY -f x'y — yX — y'x = xa + ya, (2.6) 


Ricci equations 

60 ( 0 ;) + efla) = Xa + xa' -|- Fa -|- ya' = E — |(/x -|- 3p) + (2.7) 


eo{(3) = Xfl + xj3' = -\E - |(/i -k 3p) + an - , (2.8) 

eflfl) =Yi5+ y(3'= k{( 3 - a), (2.9) 

eo{K) = Xk + xk' = (a — k)/?, (2-10) 

ei(K) = Yk + yC = —\E + — afl, (2-11) 
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Bianchi equations 

eo(-E + ^/i) = X + x(^E' + ^/i ) = — (3E + /i + p)P, (2.12) 

ei(ii/ + ^/i) = Y ^E + + ?/(£'+^/i)=3£'K. (2.13) 

Here we have used the notation / := dtf and /' := drf. 

Observe that the cyclic equations R^[ijk] = 0 are already satisfied by our choice of 
Riemann tensor. 

It might be preferable to replace the electric part E of the Weyl tensor by the 


Gaussian curvature K:=E + p,/3 + tC — of the orbits { dt = dr = 0}. The Ricci 
identities then become 

69 ( 0 ;) + ei(a) = K — ^(yU + p) — + /3^ , (2.14) 

eo(/9) = + an — |/9^ , (2-15) 

eflfl) = k{( 3 - a), (2.16) 

edn) = {a- k)(3 , (2.17) 

6 i(^) = —+ i/i + — i/3(2a + /3) . (2.18) 

The Gaussian curvature K satisfies the equations 

eo{K) = Xk+ xK'=-2I3K , (2.19) 

eflK) =Yk+ yK'= 2 kK , (2.20) 

which can be used to replace the Bianchi equations. 


3. General considerations 


In this section we discuss the choice of coordinates, the algorithm for solving the 
equations in the generic case and give an example of where essential information about 
a geometry is hidden in the coordinate gradients. 

First maybe a word about the meaning of sufficient and necessary conditions is at 
hand. The equations (|^) - (|2.13|) are sufficient, i.e., a solution to them really describes a 
geometry, but the consistency equation for, e.g., (p.8|) - (p.9|) could give “new” equations 
(in this paper one of the twice contracted Bianchi identities) that are helpful in solving 
the set of equations. 

In general, there exists several ways to choose ones coordinates. Components of 
the curvature tensor, the rotation coefficients, and/or the coordinate gradients may be 
used as coordinates. 

If we choose some of the rotation coefficients as coordinates, this is in line with 
the integration procedure in the GHP formalism [2^ suggested by Held [2^, and 

developed further in several papers by Edgar [^] and Edgar and Ludwig |[3^-|3^-li The 
commutators will generate Bianchi type equations, i.e., equations containing derivatives 
of the Riemann tensor components. In general there may be new equations among these 
(see section 0). 

§ In Edgar also pointed out the necessity of the commutator equations in the GHP formalism. 
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The theory developed in now assures that a solution to the commutator equations 
the Ricci equations (^?^ - (| 2 . 11 |) (or equivalently equations (|2.14|) - (p.l8D ), 
and the Bianchi equations (|2.12D - (|2.13D (or equivalently equations (|2.19|) - (|2.20|) ) 
describes a LRS perfect fluid solution to Einstein’s equations. Since some of the 
quantities appear differentiated in more than one equation (k, (3 and K) it could be 
helpful to calculate the consistency conditions that these equations give rise to. These 


are 


[^i! ^o] / — (r^oi ~ r^o) 6fc(/) — o6o(/) + q!6i(/) 

where f is k, fl or K respectively. By acting with eo and ei on ([^ )- (|2.11| ) one then 
gets 

eo(/i) = - (a + 2/3) (/i+ p) , (3.1) 

eflp) = a{n + p) , (3.2) 

i.e., the twice contracted Bianchi equations (note that these are not contained in 
equations ( |2.19| )-( |2.20|) ), whereas the commutator relation acting on K is identically 
satisfied. 

Let us first look at the generic case when two functionally independent quantities 
can be found among k, /3 and K. As an example we choose /3 and k as the coordinates 
t and r. Equations (p.l5D - (p.l8|) then give 


X=-\K 
Y = — a) , 

X = {a — , 

y = —hK + i/i + 


ip + 1^2 _ 3^2 


j/3 (2cr + /3) . 


This choice can be made as long as Xy — Yx 7 ^ 0. The commutator equations 
( p.5|) -( p^ are now identically satisfied due to (|2.14|) - (|2.18|) and (|3.1|) - (|3.2|) (naturally, 
since the commutators on t and r now coincide with the commutators on [3 and n). 
We hence may use the twice contracted Bianchi equations ( p.l|) -( ^^ instead of the 
commutators ( p.5|) -( p^ and we get the same set of equations as in IjT^ (our set of 
equations (p.7| )-( 2.13 ) and (|3.1J)-(|3.2|)). 


Using equations (|2.15| )- (|2.18|) for defining X, x, Y and y there now remains five 
differential equations to solve, (|2.14| ), ( p.l9|) -( p.20| ), and (pTD-(pT^. Erom (|2.14| ) we see 
that one of a and a, say a, may be freely specified. The system is then integrable since 
the commutator acting on K is identically satisfied and we just have one equation for 
each of a, /i and p. More natural would of course be to specify an equation of state 
P = pip)- Instead of equation (|3.2|) we then get 

eflp) = a{p + p) . (3.3) 
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Figure 



The consistency condition for eqnations (^.11) and 

-1 / 7 \ -1 


(as eqnations for p) is then 


efla 


eo(a) 


dp 

dp 


= 2 /s:(q! — fl) — aa 


dp 

dp 


+ a{a + 2/9) - a{a + 2/9)(/r + p) 


dp‘^ 


dp 

dp 


-2 


(3.4) 


Equations (p.l4|) and (p.4|) may be considered as developing equations along Cq for 
a and a respectively. They may then be specihed on a three dimensional space-like 
hypersurface (since a and a are specihed on this surface, their gradients along ei, which 
is tangent to the surface, may be calculated) (cf. e.g. j^). We hence have the following 
problem to solve (see hgure 1): Given an equation of state p = p{p), the development 
of p, K, a and a is given by 0 and (1^) , (|2.19|) and (| 2 . 20 |) , and (|2.14D and (|^) 
respectively, with a and a specihed on a space-like hypersurface orthogonal to Cq and 
p and K specihed on a two-dimensional orbit contained in the hypersurface. Note that 
all these considerations are purely local, since in a general space-time a foliation with 
space-like hypersurfaces may not be possible. 


3.2. A geometry that is not completely determined by the the Riemann tensor and the 
rotation coefficients 

In the example above the commutators (p.5|) - (p.6|) could be neglected when the twice 
contracted Bianchi equations where added, but cases when this is not possible do appear. 
Assume that all the Ricci rotation coefhcients and all the components of the Riemann 
tensor only depend on one functionally independent quantity (which is chosen as the hrst 
coordinate). Nevertheless the geometry could be truly depending on two coordinates, the 
second one appearing hrst in the hrst covariant derivative of the Riemann tensor. It must 
then be found in one or both of the vectors Cq and ei, since the hrst derivative should 
be constructible from the vectors, the Riemann tensor and the rotation coefhcients. We 
here give an example of such a metric. Assume that K = k = a = 0, that /9 and 
p = —3ff^ are constants and that p is a function of a (chosen to be the hrst coordinate). 
X is then found to be X = 2/3^ — — afl. Even if we add all of the Bianchi equations, 
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we will not get any more information, and we still have no equations for Y, x or y. 


But choosing them completely freely would contradict the commutators ( |2.5|) -( |2.6|) and 
we hence cannot ignore these equations. If we choose Y as our second coordinate, 
X = —(3a + /3)F from equation (|2.5|) , and (|2.6| ) gives a differential equation for y. 

—3Y^ — (3a + (3)y — y (2/3^ — a^ — a/3) + //'(3a + (3)Y = ya , 

with solution y = {C — 3a)y^/ (2/3^ — o? — a/3), where C is a constant of integration. 
The coordinates a and Y will be linearly independent as long as Xy — xY 0, 
corresponding to C + /3 7 ^ 0. If one now calculates the set one finds the second 
essential coordinate Y in the first covariant derivative (e.g. i?oioi;i = —flY) showing 
that the metric really is depending on two essential coordinates. Furthermore, the 
constant C appears first in the second covariant derivative (e.g. = —«^/3 — 

a^/3^ + ‘iaflY'^jX — CflY’^jX) showing that metrics with different C are geometrically 
inequivalent, and that the second derivative is needed to see this. 

The 1-forms for this metric are given by (see also section : 

C - 3a , 1 






{C + /3) ( 2/32 
3a + /3 


da -|- 


a^ — a/3) 
2 / 3 ^ — of 


da 


a/3 


{C + fl)Y 
dY , 


dF 


(jj 


(jj 


{C + (3)Y (C' + /3)F2 

= \Y\-l^dC+0) _ ^|^(3/3-C)/[3(/3+C)] ^ 2^|^(6/3+C)/[3(/3+C)] ^ 

3 ^ |y|-/3/(C+/3) _ ^|^(3/3-C)/[3(/3+C)] ^ ^^^~^{6y+C)/[3(p+C)] _ 


4. Static models 

In this case we have no dependence on the coordinate t, thus we have X = F = 0. The 
functions /3, n and K then have two differential equations each, so that we may obtain 
algebraic relations for some of these. 


4 . 1 . Static models with K 0 


UK 7 ^ 0 we may solve for /3 as /3 = —xn/y from equations (|2.19|) and (| 2 . 20 |) . 
Inserting this definition into equation (p.l5|) , we get an algebraic expression for, say, 
K. Combining equations (|2.17|) and (|2.18|) and substituting the algebraic equation for 
K, we find a;(p + p) = 0, i.e., a; = 0 With a: = 0, the commutator (|2.6|) gives a = 0, 


and from the Bianchi equation (|2.19| ) we see that /3 = 0. Thus the static models with 
non-planar symmetry are expansion and shear free. From (|2.15|) we have 


K = —p Y ifl' + 2aK , 

and the remaining differential equations are 
ya' = 2ati + — |(p + 3p) , 

yn' = -an + + p) , 

yp' = a(p -h p) . 


(4.1) 

(4.2) 

(4.3) 

(4.4) 
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Since we can choose one of the functions as our variable r, we see that there remains to 
solve two coupled, ordinary differential equation. 

When we have constant pressure, equation (ff.4|) gives a = 0, and equation (|4.2|) 
gives /i + 3p = 0. The only remaining equation is ( |4.3| ). If n is constant we have a space- 
time homogeneous case, and if k is non-constant we may choose it as our coordinate, 
which implies y = — p. 

When there is a non-constant pressure, it is natural to take the pressure as our 
coordinate. From equation (|4.4|) we find y = a(/i p). The remaining equations are 

(p -|- p)aa' = aiflK + a) — + 3p) , 

(/i -|- p)aK' = — a) p) . 

Most constraints on these equations will reduce the problem to an algebraic one. For 
instance, if a = k, we find that p = p + A, where A is an integration constant, and 
tC = (3p + 2A) /3. 

We may of course use another set of equations than the one chosen above. Instead 
of eliminating K with the help of equation ( |4.1|) , we may focus on eliminating a as 


K 


a = 


p 


K 


2k 


(4.5) 


This is especially useful when treating spherically symmetric space-times, i.e., K > 0. 
We define our radial coordinate r through K := r“^, so that from equation (|T^) we 
have 


y = —rn . 


Using equation (|4.1| ) to eliminate p we find that equation ( [4.3|) becomes 




pr 


dr r 

which may be integrated in terms of r and p: 

o, , dvrr — m 
K ir) = 


dvrr^ 


where we have introduced the mass function 

nr 

m{r) := dvr / /r(f)f^df , 


(4.6) 


and assumed that there are no point masses at the center (r = 0) of the model. 
Eliminating a from equation by using equation (|4.5|) , we obtain 

dp {p p) (dvrpr^ -|- m) 
dr 


2r(47rr 


m] 


(4.7) 


which is the Tolman-Oppenheimer-Volkoff equation [^]. Hence, given an equation of 
state p = p{p), one needs to solve (|4.7| ) together with the definition (|4.6| ) to get the 
full solution. In the usual approach with a metric ansatz there will remain one further 
equation to be integrated for obtaining the full line-element. Here we see that this is 
not necessary for a complete local description of the space-time. 
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We may also take another route in analyzing the equations for static LRS class 
II space-times. Since we only have dependence on one variable, and we only have one 
coordinate gradient, we may view y as a lapse function. Thus by defining dr := ydR, 
equations (f4.2| )- (^I.4| ) becomes an autonomous system of equations: 

da 21 / . dK 21 / \ dp 

— = 2aK + a - -(/i + 3p) , — = -an + k + 2 (/^ + p) , ^ = a(/^ + p) • 

This set of equations is the starting point for a dynamical systems analysis, although a 
change of variables is in place to acquire dimensionless quantities and a compact phase 
space. 


4 . 2 . Static models with K = 0 


When K = 0, the Bianchi equations can no longer be used to reduce the set of variables 
and equations. Instead, we will have 


yx' — xy' = xa + ya (4.8) 

xa' + ya' = + p) — + a'^ — a'^ + (5“^ , (4.9) 

X/d'=-|p + + OK - |/3^ , (4-10) 

yfl'= k{(3 - a) , (4.11) 

xk' = {a — K)fl , (4.12) 

yC = — afl — . (4.13) 


Now, if X = 0, the commutator equation 
gives (i) a = K, or (ii) fl = 0. Case (i) gives y 
with a cosmological constant. Case (ii) results 

ya' = —4/i — |/s;2 — an + a'^ , 
yC = , 


(^] 8 |) becomes a = 0. Equation (|4.12| ) 
-|- p = 0 , which is equivalent to vacuum 


m 


(4.14) 

(4,16) 


and 


p = + 2an . 

If we choose k = r ^ y = pi/2 + 3tfl/2 and equation ( [4.14| ) becomes an equation for pi, 
while a is an arbitrary function. 

If X 7 ^ 0, i.e., the metric is manifestly non-diagonal, we can replace equations (|4.10|) 
and (|4.13|) with two algebraic equations for p and /i: 


2,0 0/32 2 x { i 3 - a)K 


p = K + 2aK — 3/9 — 


pi = —3k^ -|- f3{2a + fl) + 


y 

2y{a — K,)fl 


The remaining equations are 

yfl' = i/3- a)K , 
xC = (a — k,)/3 , 


X 
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/ , / y{a- - a)K 

xa + ya = -1- 


X y 

,2 ^,2 I 0/32 


-an + + 2(3 — a(3 , 


X y — y X = xa + ya . 


There are in total 6 functions. Choosing one of the functions as coordinate, we determine 
X 01 y. Then there are 4 functions and 3 equations left, thus one of the functions is 
arbitrary (as in the case a; = 0). 


5. Spatially homogeneous models 

Here we have no dependence on the coordinate r. This will, as in the static case, give 
some algebraic equations from which we may express suitable quantities. 


5.1. Spatially homogeneous models with K ^ 0 

When 7 ^ 0, we use the Bianchi equations (|2.19|) and (|2.20|) to obtain k = —Y(3/X. 
From equation ( |2.18|) we then obtain an algebraic equation for K. Using equations 
( p.l5| ) and (|2.16|) is is straightforward to show that Y{y + p) = 0, i.e., U = 0. 

The curvature of the 2-surfaces {t = constant, r = constant} is related to the scalar 
curvature of the 3-surfaces of homogeneity through the Gauss-Codazzi equation. 
Using this, we find that 

(3)/? = 2K . 


If U = 0, the commutator (|2.5| ) gives a = 0, and from the Bianchi equation (|2.20|) 
we get K = 0, while ( |2.18|) becomes 

fi = K + 2a(3 + (3“^ . (5.1) 

The remaining equations are 

Xa = — a(3 — a"^ + , (5.2) 


Xfl = -\K - 1/3^ 
Xk = -2(3K . 


\P 


(5.3) 

(5.4) 


It turns out this system of equations can be reduced to one Riccati equation. First we 
note that if iF is a non-zero constant, equations e-o gives /3 = 0, p = —K ^ 


p. -|- p = 0 by ( 1^) . If iF 7 ^ 0, we can choose K as coordinate (so that X = —2/9iF by 
( ^.4|) ), and equation (p.3|) becomes linear in (3"^. Further, if we assume that p = p(iF) is 


a given function, this linear equation may be integrated: 


<-K 


kiK) = -K+ / p(iF)iF-®/2 diF . 


Thus the equation for a becomes 
da 


dK 


= F{K)G{K) + l3G{K)a + G{K)a^ , 


(5.6) 
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where F{K) := {p — (3“^ — K)/2 and G{K) := 1/{2(3K). The Riccati equation (|575D can 
be solved exactly when, e.g., there is some algebraic relationship between the coefficients 
F and G. Such a relation will in general lead to an integral equation for the function p in 
terms of K. Thus, although the choice of variables leads us to a single equation for these 
space-times, the structure of the remaining equation leads to difficulties concerning the 
equation of state. 

A perhaps physically more suitable choice of variables is 0 (expansion), a (shear), 
and p (where we have used that p and K are interchangeable, by the algebraic equation 
( ^T1) ), which leads to 

X0 =-i02-bcT^ - i(/i + 3p) , (5.6) 

Xa =+ 3a^-3Qa +p , (5.7) 

Xp = -Q{p + p). (5.8) 


Assuming non-constant density, we can solve equations (|5.6| ) and ( b-7|) for 0 and a, 
respectively, in terms of p, if we have a barotropic equation of state. 

We can look at the interesting case when the shear is small, i.e., when the deviations 
from isotropy are small. This translates into cr/0 -C 1. Linearizing the equations with 
respect to this parameter, we can solve for the expansion and shear in terms of the 
density: 

02 = F^^^p) ^ 


j,-2/3(^) A + 3p(A) j-+^ 


a = F{p) 


F-\p)- 


p + p{p) 
Q^{p) — 3p 


(5.9) 


dp + B 


9Q{p){p + p{p)) 

where A and B are constants, and F{p) := exp [J^{p + p{p))~^ dfi\. 


5.2. Friedmann-Lemaitre-Robertson-Walker models 

These space-times are characterized by the existence of a 3-dimensional spatial isotropy 
group and a 6-dimensional isometry group acting multiply transitively on 3-dimensional 
orbits. Therefore E = a = a = Y = 0, and k is a redundant variable. (Since k is no 
longer part of the set S', it is possible for it to depend on some spatial variable). This 
of course means that there is need for additional Bianchi equations, and we need fewer 
Ricci equations. It turns out that the additional equations are just the twice contracted 
Bianchi equation (|3.1J ). The redundant Ricci equations are the ones containing k. Thus, 
( ^.2|) and (|5.4|) are the necessary and sufficient equations. Choosing p as coordinate and 
an equation of state p = p{p), the system is easily integrated to give equation ( |0|) . 


5.3. Homogeneous solutions with K = 9 

When iL = 0 it is possible to have tilted models, i.e., models where the matter flow is 
not orthogonal to the surfaces of homogeneity. This means that Y ^ 0. 

The situation is similar to the plane symmetric static case. In fact, the set of 
equations becomes the same as equations (|4^) - (|4.13|) , with the replacements x,y ^ 















Construction of LRS perfect fluids 


13 


X,Y and {8/dr) (d/dt). For the generic case with Y ^ 0 the treatment is identical 
to the one in the static case. For Y = 0, one obtains a = 0 and — a) = 0, i.e., 
K = 0 or a = fl. In both cases the systems can be completely integrated, similarly to 
the static case. 

If one assumes dust, i.e., p = 0, the equations may be completely integrated as was 
shown already in [^]. If n is used as time coordinate, one gets 

X = =FK^\/l + Ck , Y = tC -^ , fl = ±K\/l + Ck , 

3Ck 

^ p{3Ck + 2) ^ _ 3CDk^ 

^2VTTV^ GCk^VTTC^ ’ d-kVTTc1^ ’ 

where C and D are constants of integration. There is also an isolated solution p = 3CrC, 
but then y = 0. If k cannot be used as coordinate, i.e., if it equals a constant, it also 
follows that y = 0 . 


6. Inhomogeneous, time dependent models 


Herein lies perhaps the most interesting models, such as collapsing perfect fluids with 
an inhomogeneous matter distribution. Since no special simplifications appear in this 
class, the system of equations is in general impossible to solve explicitly. One therefore 
has to resort to simplifying assumptions. Some special cases are treated below. 

No particular simplification is implied by the constraint iF = 0, since the Bianchi 
equations are identically satisfied by this choice. The analysis of the system of equations 
characterizing these space-times can be performed analogously to the case iF 7 ^ 0 . 

Before turning to the examples we show how our set of equations may be written in 
a form similar to that one would obtain from a metric ansatz. In the generic case with 
TF 7 ^ 0, we may proceed as follows: Define a and a through equations ( p.5| ) and (|2.6| ), 
thus expressing them as derivatives of the coordinate gradients. Further, define f3 and k 
by equations (|2.19| ) and (|2.20| ), respectively. From equations ( |2.15| ) and (|2.18| ) we And p 
and p respectively. We are then left with the second order equations ( p.l4| ), ( p.l 6 | ), and 
( P.17| ) for the variables K, x, y, X, and Y. It is straightforward to show that equations 
( p.lfl ) and (I2.17D are the same. By performing a coordinate transformation t' = t'{t,r), 
r' = r'(t, r) it is always possible to get y = 0 and a; = 0 , hence making d/dt' parallel to 
the four-velocity of the fluid. The remaining freedom in choice of coordinates are now 
transformations of the type r" = r"(r'), t" = t"(t'). With this choice of coordinates, 
equations (p.l 6 |) and ( p.l4| ) reduce to (dropping all the primes) 
d'^f dx df dv df 1 df df 


dx df dv df 

dtdr ^ dr dt dt dr 


2 dt dr 


2 e« - 


2v 


—e 


2 X 


d^f d^X 

+ 2 — - 2 

dr-i ^ dr^ 

d^f d^v 

— - 2 -h 2 

dC dC 


dr 

dv 

'di 


= 0 


dx d^ dv df 2 

dr dr dr dr dr dr 


( 6 . 1 ) 


dt 


dt dt dt dt dt dt 


= 0 , ( 6 . 2 ) 
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and 
, If 


respectively. Here f := \nK, x •= InX and v := Inr/. 

These are essentially the equations one would obtain from a metric ansatz, 
since the functions X, y and K correspond to the metric coefficients (see section 
Non-comoving coordinates has been used by some authors (see ||l9|-[^ 
references therein) to analyze spherically symmetric perfect fluids (see also 
the line-element is to remain diagonal, we must have xX — yY = 0 (see section |^. 
In general, this is a too weak constraint to yield any significant simplifications of the 
equations, and further assumptions has to be made. 

We now continue presenting some of the special cases appearing in this class of 
space-times. 

The shear- and expansion free fluid space-times have a = fl = 0, and they 
are therefore a bit pathological because of their time-dependence (although the static 
subcase is of great physical interest). In fact, the equations ( |2.14|) -( |2.18|) and ( p.l9| )- 
( p.20| ) reduce to 

efla) = \K — |/r — — aK-\- of , (6.3) 

efln) = -\K + \pL+ , (6.4) 

eflK) = 2 kK , (6.5) 

(which are just the static set of equations when x is zero), with p = —K -|- 2aK tfl, and 

eo(p) = 2 ^ 60 (a) , 

eo(«) = eo(X) = eo(/i) = 0 , 

where the last equality follows from (|3.1|). If we align our spatial coordinate so that 
a; = 0, we see that all the time-dependence appears in p, a, and the gradients of the 
coordinates. This is a manifestation of the fact that this type of solution may be 


generated from a static solution ||26 


First, suppose that K' = Q. Then it is straightforward to show that either iF = 0 
or /i -|- p = 0. Therefore, if X > 0 (the case iF < 0 can be treated in an analogous 
manner), we choose K := 1/r^ x = 0, p = —tk, and we And that the commutator 
( ^^ is identically satisfied. Assuming that p is a given function of r, we may then 
integrate equation (1^ ) for k. Alternatively we may define p through equation ( |6.4|) . 

If eo(a) = Xd 7 ^ 0, we can choose a := f, and Y is given by equation (|6.3| ). Equation 


( p.5|) can then be integrated for X, since all other functions are known. 

If Co (a) = 0 (-v^ eo(p) = 0), all the rotation coefficients and curvature components 
depend solely on r. Thus all the time-dependence enters through the coordinate 
gradients. This is another example of when the coordinate gradients contain essential 
information about space-time. If X 7 ^ 0, we take t = X, and Y can be obtained from 
( p.5| ) since a is in principle known from equation ( |6.3|) . If X = 0 but Y 0, the 
analogous result holds. 

Shear free fluids are characterized hy a = fl. If we transform our coordinates 
such that Y = 0, equation (|2.16|) gives j3 = j3(t). If we make one further coordinate 
transformation, to obtain x = 0 , we can use equations ([^ ) and ( | 6 . 2 D as our held 
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equations. The shear-free condition can then be stated as v = f/2 + g{r), where g is an 
integration fnnction. Fnrther, eqnation (|6.1| ) can be integrated so that x = — + f{t), 

where / in an integration function. Thus, equation (|6.2|) becomes an equation solely 
for and all the time-derivatives in the last bracket drop ont. We define onr spatial 
coordinate r by choosing a special form of the function g. 


6.1. Acceleration free fluids (including dust) 

For constant pressure the equations may be completely integrated, as is indicated in 
If one chooses Y = a; = 0 it follows from (|3.2|) that the pressure is a function of t only. 


From ( p.6| ) one gets that also X only depends on t. We then still have the freedom to 
make X = 1 by performing a coordinate transformation t' = t'{f). This is in accordance 
with the possibility to choose a synchronons and comoving frame for non-accelerated 
matter ||^. Note that this choice conld contradict an arbitrary choice of coordinates 
from . From equations (|2.19|) and ( |2.17|) it then easily follows that 

K = K^fir), 

where /(r) is an arbitrary fnnction of r. The snbseqnent solving of the eqnations is 
greatly simplified if we introduce the function S by 

1 


^:= - 


«/' 


As will be seen in the last section, 8 “^ will then appear in the metric in front of dfl^. 
The remaining qnantities may then be expressed in terms of S as 

1 


y = 


fS' 

p = 2 


S' s 
S' ' ^ S' 


f 1 
pss' ^ 52 


1 SS' 52 

-h 2-^- 

7252 SS' A2 


and S is given by the following equation 

1 1 PS^ ^ 

SS+Y + 5-2P + —-0 


(see also ||2^). The pressnre p = p{f) is free to specify. If p is a constant one may solve 
for t as a function of S and r. The solution is given by 


t{S,r) = ± 



S 


dA + fi , 


S(/-" - 1) + 9 - IpS^ 
where g and h are arbitrary fnnctions of r. If p = 0 the integral may in fact be expressed 
in terms of elementary functions. Even if one cannot solve for S = S{t) in general, one 
can get an explicit line-element by changing coordinate from t to S. Off-diagonal terms 
will then of course be introduced in the line-element. 
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Kitamura found a class of solutions with non-zero shear, expansion, and acceleration, 
generalizing the results of Gutman and Bespal’ko (a solution rediscovered by 
Wesson H). Lake |@, Van den Bergh and Wils ||4^, Collins and Lang iQ, Coley 
and Tupper [i^, and Sussman @|. The solution of Kitamura can be obtained using 
the eqnations presented here with the assnmptions x = Y = 0, X = —y, a = n, and 
a = a{t). 

Using non-comoving coordinates, Knutsen [^] found a spherically symmetric 
perfect fluid solution with shear, expansion and acceleration. This solution is from 
a physical point of view quite well behaved, but as shown in |^] it can not be matched 
to a vacuum exterior. Also, it turns out that in general it is not possible to nse comoving 
coordinates to express the solution in terms of “well known analytic functions” 


24 


Another solution with all kinematic quantities non-zero, which is new to the authors 
knowledge, may be found with the ansatz 

Y = x = 0, K = f{t)/r\ 
a = 0 , a = K . 

The form of K defines the coordinate r, while retaining some freedom in the choice of 
time-coordinate. 

From equation ( |2.17| ) we see that n = ^(r), and from equations ( |2.20| ) and ( |2.5| ) we 

find 

9{t) 


y = —nr 


X = 


where g is an integration function. Further, equation (p.l6|) can be integrated so that 

r 

and inserting this into the Bianchi equation ( |2.19|) we have 

^/ + 2F/ = 0, (6.6) 

which we use to define g. Expression for p and fi can be obtained from equation (|2.14|) 
and (|2.15|) respectively: 


3/ . 

P = ^ - 2/i 


+ 


3F2 


2 / 9 

11 = — - 3 k ^ -k 


3F2 gF 
——I-:r 




Inserting this expression for y in equation (|2.18|) , we can integrate and find 

where U, B and C are constants. With this, y = p + QB. f may be chosen as the time 
coordinate t. 

We also see that y = Cp/r'^ + 3B, so that y is monotonically decreasing with 
increasing r, as is p. The fluid as such is thus reasonably well behaved, although 
it harbours a singularity at r = 0, but unfortunately it can not be matched to the 
exterior Schwarzschild solution. The reason for this is the following: Among the junction 
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condition (using the Darmois-Israel matching procedure ||47|,|4^) on the matching surface 
are a = Oext and k = /text, where we have used the subscript “ext” to denote the 
corresponding quantities for the exterior solution. For the particular fluid solution 
presented here (and also for the solution of Kitamura), we have a = k. Since this 
should hold on the junction surface, we have Oext = «^ext on that surface. If we denote 
Scwarzschild’s radial coordinate by R, the latter condition gives a constant matching 
radius R = 3m (where m is the mass parameter of the Schwarzschild solution). On the 
other hand, from the remaining matching conditions on flnds that the matching radius 
is a time-dependent quantity, thus giving a contradiction. This problem also arises in 
Knutsen’s solution. 

It is straightforward to show that Kitamura’s solution and the solution presented in 
this section are not generically equivalent. Equivalence can be obtained for certain values 
of the integration constants, but with arbitrary values of the integration constants in the 
two solutions there is no equivalence. Also, our solution is not equivalent to Knutsen’s. 


6.3. Space-times with a homothetic motion 


If space-time is self-similar, i.e., has a homothetic Killing vector f the held equations (for 
the LRS case) are reduced to ordinary differential equations (see |^] for a review). This 
is not the only reason for studying these space-times. They often contain interesting 
physical phenomena like, e.g., shock waves, black hole formation, evolution of voids, 
possible violation of cosmic censorship, etc. mM- Space-time has a homothetic Killing 
vector if 

where the constant c may be chosen to c = 1. If the coordinate u is adapted to the 
homothetic vector, so that ^ it is easy to show that all components of the metric 

will depend on u as ( 7 ^,^ = where ( 7 °^ is independent of m. It then follows that the 

components of the frame vectors Cj depend on u as e““. The Ricci rotation coefficients 
■jijk also depend on u as and the Riemann tensor Rijki as (see also ||58| ). 
Inserting this into the IC (|2.5| )- (|2.13| ) gives a system of ordinary differential equations. 
Note that we no longer can put x = Y = 0 in general, since this could contradict 
the choice f = du. However, it is still possible to perform a coordinate transformation 
u' = u f{v), where v is the other essential coordinate, giving f = du'. Hence the 
functional dependence on u' will be the same as on u. With such a transformation it is 
possible to remove Y in the time-like case and x in the space-like. 

Concerning the equation of state we have the following well known result: If we 
assume that the equation of state is barotropic, i.e., p = p{p) (which is reasonable in 
the case of a perfect fluid), the homothetic motion gives p = (7 — 1 ) 71 , where 7 is a 
constant. 
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Temporal self-similarity Here t = u. In general the system still gets quite complicated: 

y 

! r? aj3 

K = -^-, 

y X 


xa XK, ^ a j3 

y2 y2 y y 


,y 


xa' + ya' = y — S—afl -\- — ok — 2af3 , 


X 


supplemented by the defining equations 

X=fl+—, 

y 


.y 


K = jj, -\- — flifla + fl) — 2 — 0/9 , 

X 


7 


y 


X 


/i = — I OK + afl H— a(3 H —na 


X 


y 


7 


+ a-y = 2k- 2-(7 - l)/3 . 

7 — 1 a; X 

When rewritten in terms of dimensionless variables (chosen such that the phase 
space is compact), this system (and the corresponding spatial one) has been studied as 
a dynamical system in |^,|^ (see also 0)- 

We here look at some simple special cases. The assumption a = fl (no shear) gives 
Po := pC = constant, a = 0, W = and p = S/dge”^*. If one chooses the coordinate r 
as r = x/y one gets 


K = -^ = = (37 - 2)^0 

rj^ rjrt‘2^ 2 ?^ 


xCJ2slL,->, 


The pressure is given by p = (7 — 1)/U and all values in the interval 1 < 7 < 2 are 
allowed. 

The simplification a; = 0 leads to W = 0 (otherwise one gets /i = 0). If one chooses 
the space-like coordinate r as r = pC the solution is given by 


X = 


27 


2-7 


-re 


— t 


a = 


2-37 

y = 

—t 


a = 


27 

-re 


37 -t 


7 


KQ = r 


2-7 

2 Dr^ , 


P = re 


-t 


and 


p = (7 - 1 )/U = ^ — -Dr^-^e , 

2-7 


where D and 7 are constants. We see that the allowed interval for 7 is given by 1 < 7 < 2. 
The 1-forms are given by (see also section on how to obtain and 




2 — 7 e* 1 2 — 70* , 

dt , -dr , 


27 r 


2 — 87 Kgr 


7-2 2-7 ^ 7-2 2-7 ^ 

uj = r 2-37 6 27 dr , uj = r 2-370 27 dtr . 
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Table 1 . Definitions of the quantities used in section 

Case a> F k ^ 

Static a;°| /3 It, time-like 

Spatially homogeneous —k 1 r, space-like 

Inhomogeneous and time dependent — 0 0 — 


Note that in these coordinates all the components g^y of the metric are not proportional 
to This proportionality may however be achieved by a coordinate transformation of 
the non-essential coordinates v and w. 


Spatial self-similarity Here r = u. These space-times can be analyzed analogonsly to 
the temporal case. As an example of an exact solution, we look at the case when T = 0. 
This implies k = a = —y, a = 0, and p = p. We find that k is a constant, and 
K = p + tC — . If /ie^^ is chosen as the time-like coordinate one obtains 

= jjL + e~^ , 

where C is an integration constant. 


7. The isometry algebra 

The isometry algebra may be obtained by projecting the 1-forms onto the orbits by 
putting x°‘ = constant and all generators of the Lorentz group, rf, except t% to zero. 
Defining := uj% = t\ -|- 7 ^ 3 * 0 ;®, the isometry algebra on the orbits is given by (by | 
we denote projection onto the orbits) 


k dcj 

= 0 , 

(7.1) 

du:^\ 

= a -I- F(jj A , 

(7.2) 

dcj^l 

= u>^\ A -|- Fu) A , 

(7.-3) 

doj^l 

= Ku;^\Auj^\ , 

(7.4) 


where lj, F, and k are defined in Table 1. 


7.1. Spherical and hyperbolic symmetric space-times {K 7 ^ 0) 


When K ^ 0 the integrability conditions give F = 0 (see sections | and |5|) . Depending 
on the value of sgn(F) the algebra (|7.2|) - (|7.4|) has two different solutions: 


(a) sgn(F) = 1. 

(— sin C d^ -I- sin 6 cos ( d0) , 
I = (cos C -I- sin 6 sin ( dcf) , 

= cos 6 d0 -1- dC • 
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(b) sgn{K) = -1. 


uj'^l = \K\ sinC d6'+ sinh6'cosC d0) , 

(cosC d^ + sinh6^ sinC dcf) , 


LJ 


= cosh 6* d0 + dC . 


For the static and SH cases, equation (|7.1|) has the solution u) 
The line elements are then given by 


= de 


Static case 


where 


9 = 


dC - y-^ dr2 - iF-i ( dO^ + sin^ 9 dcf^) , 
dt^ — y~‘^ dr^ — ( d6^ + sinh^ 6 dcffl') 


A{r) = , 


with Aq a constant; 


SH case 


where 


9 = 


X-2 dt^ - dr2 
X -2 - 52 dr 2 


B{t) = , 


with Bq a constant; 


iF-i (d02 + sin2 0d02) , 
\K\-^ {de^ + sinlfl 9 d(j)^) 


Generic case 

A-2 [(|/2 - x^) dfl + 2(a:X - yY) dt dr - {X^ 
_l -K-^ {d9^+ siifl9d(j)^) , 

^ A“2 [(y 2 _ ^ 2 ^ _|_ 2(^xX — yY) dt dr — (A^ 

-|A|-i (d02 + sinh2 0d02) , 

where A := yX — xY. 


K>0 
K <0 


K>0 
K <0 


y2) ^^2] 
y2) ^^2] 


A > 0 
A < 0 


7.2. Plane symmetric space-times {K = 0) 

If a; = d.^, then the algebra (|7.2|) - (|7.4| ) has the solution 
ufl\ = (cos C dr — sin ( dtr) , 

(sin(C dr + cosC dtr) , 

= dC , 

where r and w are space-like. Putting F = 0 we obtain the solution to the generic 
algebra. 

The metrics become 
Static case 

9 = A^ dt^ — y~‘^{ dr — xA dt)^ — ( dr^ -|- dtr^) , 
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A{r) := , 

/ r 

- df + Dst , 

y 

with Cs and Dg constants; 

SH case 


g = X-\ dt - YG drf - dr^ - ( dn^ + dw^) , 

where 

G{t) := , 

H{t, r) := - j ^ di + DhV , 


with Gh and Dh constants; 

Generic case 

g = A“^ [(r/^ — x^) dt^ + 2(a;X — yY) dt dr — (X^ — dr^] 

— (dr^ + dtr^) , 

where / = f(t,r) satisfies 

eo(/)=/5, ei{f) = -K, 

which is a set of integrable eqnations dne to the integrability conditions. 
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